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ABSTRACT 

Extinction in galaxies affects their observed properties. In scenarios describ- 
ing the distribution of dust and stars in individual disk galaxies, the amplitude of 
the extinction can be modulated by the inclination of the galaxies. In this work 
we investigate the inclination dependency in composite spectra of star-forming 
disk galaxies from the Sloan Digital Sky Survey Data Release 5. In a volume- 
limited sample within a redshift range 0.065 — 0.075 and a r-band Petrosian 
absolute magnitude range -19.5 to -22 mag which exhibits a flat distribution of 
inclination, the inclined relative to face-on extinction in the stellar continuum is 
found empirically to increase with inclination in the g, r and % bands. Within the 
central 0.5 intrinsic half-light radius of the galaxies, the g-band relative extinction 
in the stellar continuum for the highly-inclined objects (axis ratio b/a = 0.1) is 
1.2 mag, agreeing with previous studies. The extinction curve of the disk galaxies 
is given in the restframe wavelengths 3700 — 8000 A, identified with major optical 
emission and absorption lines in diagnostics. The Balmer decrement, Ha/H/3, re- 
mains constant with inclination, suggesting a different kind of dust configuration 
and/or reddening mechanism in the H II region from that in the stellar contin- 
uum. One factor is shown to be the presence of spatially non-uniform interstellar 
extinction, presumably caused by clumped dust in the vicinity of the H II region. 

Subject headings: galaxies: fundamental parameters — techniques: spectroscopic 
— methods: data analysis 
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Introduction 



Extinction in galaxies affects their observed properties. Ideally it should be taken into 
account when inferring the intrinsic properties of the galaxies. One empirical approach in 
studying extinction in disk galaxies is to compar e the observed p roperties of objects with 
similar intrinsic properties at different inclination ( Holmberg 19581 ). Optically thin galaxies, 



when inclined, should show higher surface-brightness than their face-on counterparts because 
of a larger column of stars observed in a smaller projected area on the sky. This approach 
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Recently, there are revived interests on the topic of extinction because of the large 
sample of gala xies made availa ble by various sky surveys including the Sloan Digital Sky 
Survey (SDSS; I York et al.l 120001 ) . Progress has been made in studying extinction in galaxies 
empiricall y by using t e chniq ues that are depend on or independent of the inclination-related 
approach. IShao et al.l ( 120071 ) demonstrated the inclination dependency of dust extinction in 

and an inclination 



disk-dominated galaxies from the SDSS, by using the z- band extinction^ 



measure being the apparent major to minor axes ratio. iHolwerda et al. 



(120071 ) applied the 



technique of occulting galaxy pairs ( [White fc Keellll992h to 83 SDSS spiral disk galaxies in 
in probing the opacity of the foreground galaxies, and found that the optical depths (r) in 
the g, r, i bands are typically within the range — 1.5 across the di s k from approximately 
0.5 to 4 half-light radii (r e ff, in their Figure 6). lUnterborn &: Rydenl ( 120081 ) looked at disk- 
dominated galaxies from the SDSS and found that the highly inclined objects (b/a ~ 0.1) 
exhibit ~ 1 mag extinction correction relative to the face-on magnitude. iMaller et al.l (120091 ) 



studied the observed photometric properties of the SDSS galaxies as a function of inclination, 
and concluded that the inclination-dependent extinction correction in the g band could reach 
~_lmag with a median value of 0.3 mag for a sample of disk galaxies. iBailin fc Harris 
(200 81) applied the inclination-corrected magnitudes to the context of galaxy classification. 
Driver et al.l ( 120071 ) developed a sophisticated approach to derive iteratively the extinction- 
inclination relation for both the bulge and disk components of disk galaxies, and found that 
the edge-on to face-on relative extinction is ~ 2 mag and ~ 1 mag, respectively for the bulge 
and the disk. This is an unique approach to date to simultaneously perform the bulge-to-disk 
decomposition and derive the extinction-inclination relation for both components. 



1 The z-band extinction was estimated by iKauffmann et al.l (|2003l) who used the difference between the 
colors from the stellar population model and from the measurement, and extrapolation to the z-band as- 
suming a standard extinction curve. 
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Inspired by these studies, and that the SDSS provides for us a large sample of galaxy 
spectra within a fixed central 3"-diameter fiber, and th at the calibration of the spectra is 



i ndepe n dent of that in the SP S S pho tomet ric data used b ylShao et al.l (120071 ) , iHolwerda et al. 



( 120071 ) . lUnterborn fc Rvdenl (120081 ) and Mailer et al.l ( 120091 ). we now investigate if there 
is any inclination dependency in the observed galaxy spectra; and what is its wavelength 
dependency, e.g., the effect in the gaseous component, such as the H II region, vs. that 
in the stellar continuum; and the comparison of extinctions among different elements. Our 
approach is to explore the inclination dependency of composite spectra of star-forming disk 
galaxies in a well defined volume-limited sample that exhibits a flat distribution of the galaxy 
inclination. The volume-limited nature of the sample eliminates the Malmquist bia^l in the 
analysis. The sample is also selected to exhibit a flat distribution of the inclination angle, 
so that at a chosen redshift (z) range, a galaxy of luminosity in a defined range would be 
observed by the survey, independent of the inclination angle (or other properties) of the 
galaxy. It is a very important step because one can then assume randomness in the intrinsic 
properties of the sampled galaxies. We further confine ourselves to a sample of a narrow 
redshift range, so that any difference detected in the observed properties between the two 
composite spectra at different inclination would not be mainly due to the difference in the 
redshifts of the galaxies. On the other hand, the composite spectra, or the sample-average 
spectra, are used in an effort to average the variation in the intrinsic properties of the 
galaxies, as the actual optical depth of individual galaxies depends on detailed properties 
such as the distributi on of dust, gas, and stars, the structu re of the galaxies such as the 
presence of arms (e.g. lElmegreenl Il980l ; iBeckman et al.lll996f ). the Hubble type, etc.. 



We introduce the galaxy sample in S|2j We present results on the inclination dependency 
of composite spectra of the galaxies in §31 We derive the empirical continuum extinction in 
the SDSS g, r, i bands in §U the model-based face-on extinction in §4.31 and the extinction 
curve (the extinction as a function of restframe wavelength) in §4.41 We study the relation 
between the Balmer decrement and the inclination- dependent continuum extinction in §5J 
We summarize the results in §6j 

We use the term "reddening" to refer to the reduction in flux in the shorter wavelengths 
relative to longer wavelen gths, or a decrease in the steepness of the spectral slope. We follow 
Davies &: Bursteinl (119951 ) and use the terminology "extinction" to refer to the reduction of 
flux due to absorption and/or scattering by dust. We use "attenuation" interchangeably 
with extinction. Following the convention in the SDSS we express the spectra in vacuum 
wavelengths. 



2 In that, one preferentially draws galaxies of higher (or lower) luminosity at higher (or lower) redshift, 
for example, in a flux-limited survey like the SDSS. 
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2. The Sample 



As part of the Sloan Digital Sky Survey (jYork et al.ll2000l ) spectra are taken with fibers 



of 3" diameter (corresponding to 0.18 mm at the focal plane for the 2.5 m, f/5 telescope). 
All sources are sel ected from an initial imaging survey using the SPSS camera describ ed in 



Gunn et al.l ( 119981 ) with the filter respons e curves as described in iFukugita et al.l ( 119961 ) , and 



using the imag ing process i ng pi peline of iLupton et al.l ( 1200 ll ). The astrometric calibration 



is aescriDea ] 
(Smith et al. 
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2002; 
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galaxies for spectroscopic follow-up observations are discussed by IStrauss et al.l (120021 ) . 



Our sample is taken from the spectra of star-forming disk galaxies in the SDSS Data 
Release 5 (DR5). We determine the c r iterio n for selecting disk-dominated galaxies with the 
approach used in lUnterborn &; Rydenl (120081 ) , who characterized that galaxies with exponen- 
tial radial light profiles fall on the blue sequence in the color-magnitude diagram. Based on 
the full DR5 spectroscopic sample, we define a parent sample of galaxies at redshifts — 0.2 , 
which covers the median redshift of the Main galaxy sample, about 0.1 (IStrauss et al.ll2002l ). 
The galaxies in the parent sample have fracDeV r = — 1 in the r band, where fracDeV r equals 
for a pure exponential, and e quals 1 for a pure d e Vaucouleurs light profile as provided 
by the SDSS imaging pipeline (ILupton et al.l 1200 ll ). The color- magnitude diagram of the 
parent sample is shown in the left of Figure [TJ in the r-band Petrosian absolute magnitude, 
M r , vs. u - r. Both the absolute magnitudes and the colo rs are K-corrected, a nd corrected 



against foreground dust using the Milky Way dust map (ISchlegel et al.lll998l ). Two main 
distributions are discernible in the color-magnitude space, corresponding respectively to the 
red sequence (u - r > 2.4) and the blue sequence. Based on the parent sample, we deter- 
mine the appropriate fracDeV r cut for selecting the blue galaxies by repeatedly narrowing 
the fracDeV r range and inspecting the resultant color-magnitude diagrams. Galaxies with 
fracDeV r in the range of — 0.1 are fairly smooth in the color- magnitude distribution (the 
right of Figured]), with a contamination by the red galaxies of ~ 5% (number ratio of the 
red sequence to the total). The galaxies with fracDeV r = — 0.1 are used in the subsequent 
analysis. 

The inclination of a galaxy is represented by the ratio between the minor and major 
axes in the r band, expab r (or b/a in this paper for simplicity) provided by the SDSS through 
fitting each galaxy image with an exponential light profile, arbitrary axis ration and position 
angles. Apparent axis rati o is shown t o be a good measure of the average int rinsic inclinat ion 
of a sample of galaxies by IShao et al.l (120071 ) , who used the approach as in iRydenl (120041 ) to 
simulate the average relationship between the intrinsic and apparent inclinations in a sample 
of galaxies of Gaussian-distributed intrinsic thickness, intrinsic ellipticity that are observed 
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at different viewing ang lesE 



2.1. Volume-Limited Sample 

Based on the parent sample with a constrained fracDeV r = — 0.1, several volume- 
limited samples are constructed by selecting galaxies in narrow ranges of redshift. From these 
we seek sample(s) that exhibit a flat distribution of the galaxy inclination. We emphasize 
that this approach is purely empirical but is suitable for our purpose, as such a volume- 
limited sample with a flat distribution of galaxy inclination is obtained. By the nature of 
this approach, we neglect the step to correct for any possible biases that are likely to exist 
in the galaxy sample. Rather, we resort to the cancellation of all of the possible biases that 
gives rise to a flat distribution of inclination. In the current context, a biased sample is 
the one which preferentially includes more face-on than edge-on galaxies, or vice versa. We 
summarize three major kinds of biases: 

1. The flux-limited nature of the SDSS - the flux limit in the SDSS spectra forbids the 
inclusion of inclined galaxies that are too dim to be observed at certain redshift. 

2. The use of observed properties in sample selection - the observed properties are pre- 
sumably dust-attenuated. The use of the inclination-induced reddened colors would 
cause a bias against edge-on galaxies. B esides, the bulge- to-total flux ratio (B/T) of 



inclined galaxies was found to decrease ( fluffs et al.l I2004J ) because the bulge is more 



attenuated than the disk, implying the use of observed fracDeV r (correlated with B/T) 
would cause a bias towards edge-on galaxies. 

3. The bias that is tough to be controlled without further simulation - both the intrinsic 
thickness of galaxies, and the seeing, would make edge-on galaxies to appear rounder. 

The sample selection for studying inclination-dependent galaxy properties is therefore a 
highly complex issue. This is exactly why we seek an empirical approach as a pilot study. 
Naturally future work should be focused on eliminating these biases in the SDSS spectro- 
scopic samples that are used for inclination studies. One possible approach is to derive all 



We note that, however, the galaxy sample s they conside red are different from ours, so that IShao et al 



(|2007fl used simulated disk-only galaxies, and iRvdenl (|2004r i used disk-dominated galaxies from the SDSS 



at similar redshift range to us but larger apparent half-light radius, > 5" (see Figure [5] for the size of our 
galaxies). 
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of the observed properties of interest that are subjected to both the inclination effects and 
the selection biases, through empirical relations in the literature and/or simulation. We may 
then work backwards to derive the intrinsic properties of these galaxies. 

To continue with the sample selection we plot the number distribution of the axis ratio 
of the galaxies (b/a) in Figure 2(a) to 2(c). Solely based on redshift cuts and without 



additional magnitude cuts on top of the survey magnitude limits, we see that at lower 

0.03 — 0.04) the face-on galaxies are relatively under-sampled. 

= 0.1 — 0.11) the edge-on galaxies are relatively 



redshifts (Figure 2(a) 



Conversely, at higher redshifts (Figure 2(c 



under-sampled. The upper and lower absolute magnitude limits are illustrated in Figure |3l 
as a function of redshift. As the SDSS spectroscopic sample is flux-limited (r-band Petrosian 
apparent magnitude brighter than 17.77 mag), we deduce that the lack of edge-on galaxies at 
higher redshifts is due to the dropping out of the survey apparent magnitude limit for those 
highly inclined objects. At the redshift of 0.1, the highly inclined objects that are dropped 
out of the flux-limited survey would be, as inferred from Figure [3j intrinsically dimmer than 
~ -20. Conversely, the lack of face-on galaxies at lower redshifts is expected to be related 
to the interplay among all of the selection biases described previously, causing more edge-on 
systems to be picked. 

To understand further this result, we calculate the average inclination of the galaxies, 
plotted in the plane of M r vs. redshift (Figure H]). Only galaxies at z ~ 0.065 — 0.075 
show an average inclination of ~ 0.5, consistent with the expected average inclination of a 
random sample of galaxies that is unbiased against inclination and the range of the inclination 
being within and 1. This consistency, we note, does not prove that those galaxies are 
free of any selection bias. One such indication is that the average b/a value for brighter 
galaxies (brighter than M r fh -21) is not 0.5, but 0.65. This bias, however, does not affect 
the average luminosity in each of the composite spectra, because the average luminosity is 
heavily weighted towards the characteristic luminosity of a galaxy sample and less sensitive 
to the number of the bright-end objects. Besides, our lower absolute magnitude limit, -19.5, 
is dimmer than the characteristic magnitudes (the median M r ranged from -20.1 to -20.4 
for b/a from 0.1 to 1.0 for galaxies at z — 0.065 — 0.075 in Figure H]), meaning that the 
characteristic luminosity is included in the luminosity distribution. 

Because the z = 0.065—0.075 sample is mostly flat in terms of its inclination distribution 
(Figure 2(b) ), the relatively small number of nearly edge-on (b/a = — 0.1) and nearly face- 
on (b/a = 0.9 — 1) galaxies in Figure 2(b) is interpreted t o be physical. Discussed by other 
authors (e.g., IShao et al.l 120071 ; lUnterborn fc Rydenl 120081 ) . in an unbiased sample the lack 



of exactly or nearly face-on galaxies is due to intrinsic elli pticity of the galaxies (on average 
e = 1 — b/a = 0.16 in the SDSS disk galaxies, iRydenl 120041 ). and the lack of exactly or nearly 
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edge-on galaxies is due to intrinsic thickness of the galaxy disks. Another factor that may 
come into play is that the axis ratio b/a may not be a pure inclination measurement of the 
disk component only, but of the combined bulge+disk systems which are subjected to the 
seeing effect. This factor could explain, e.g., the lack of edge-on galaxies. In the following, 
we confine ourselves to the galaxies located at redshift 0.065 — 0.075, with a corresponding 
survey absolute magnitude limits of ~ -19.5 to -22 mag in the r band. 

Finally, only spectroscopically-classified star-forming galaxies are considered, where 
narrow-line active galactic nuclei are removed from the sample through the line ratios 
[O IIl]/H/3 vs. [N II] /Ha according to the calibration by Kewley et al.l ( 120061 ). In doing 
so the origin of the emission lines of the galaxies considered in the later session, e.g., the 
Balmer decrement (Ha/H/3), is well constrained - arise from the H II region that is excited 
by the O/B stars. We take the equivalent widths (EW's) of these emission lines as those 
measured by the S PSS, and add a constant restframe stel lar absorption of 1.3 A to each 
of the Balmer lines ( Hopkins et al.ll2003l ; iMiller et al.l 120031 ). Each of the involved emission 
lines is required to be larger than 0.7 A in the restframe EW, and is of at least 2a detection. 
These selection criteria provide higher confidence to the Balmer decrement measurements 
that are performed in §3 Composite spectra are constructed in several bins of inclination: 
b/a = 0.1 - 0.2,0.2 - 0.3,0.3 - 0.4,0.4 - 0.5,0.5 - 0.6,0.6 - 0.7,0.7 - 0.8,0.8 - 0.9, and 
0.9 — 1.0. The bin 0.0 — 1.0 is dropped in the subsequent analysis because of the limited 
number of galaxies (Table [1]). 



2.2. Fraction of Galaxy Light through Spectroscopic Fiber 



The apparent half-light radius of the star-forming di sk galaxies is found to increase with 
the galaxy inclination, as shown in Figure [5) Similar to iMaller et al.l (120091 ). the relation in 
our galaxies can be described by 



logic i r cs b/a ) = logio (^ff 1 ) - P r logio ( b / c 



where r e fj b//a is the apparent half-light radius of the inclined galaxies, and r^ 1 is the value of 
the face-on galaxies. Through a linear least-square fit, the best-fit coefficient i n the r band, 
r , is found to be 0.26 ± 0.01. This value agrees fairly well with that from IMaller et al. 
(120091 ). who found /3 r = 0.2 for a different sample of nearby SDSS disk galaxies which also 



showed inclination dependency in their observed proper ties. The larger apparent half- light 
radius in the more incli ned disk galaxies was seen also in iHuizinga fc van Albadal (119921 ) and 
Mollenhoff et al.l ( 120061 ). respectively in an independent galaxy sample and in a model for 
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galaxy spectra. The adoption of the relation in Eqn. [T] is physically motivated. Its origin is 
presumably related to the extinction radial gradient in the galaxies, so that in the inclined 
galaxies the central region dims by a greater degree than the outer regions, resulting in an 
increased apparent half-light radius. This explanation is consistent with what we found later 
(see Figured]), that dust is present at least in the disk galaxies, and more so for the more 
inclined galaxies. 

Correcting the half-light radius of each galaxy by using Eqn. [TJ the sample-average 
ratio between the radius of the SDSS spectroscopic fiber (a constant r fiber = 1.5") and the 
face-on r e fj (i.e., rgber/fee 1 ) is calculated to be approximately 0.45. This ratio stays fairly 
constant with the galaxy inclination, implying that a similar fraction of the full galaxy light 
is measured in the spectra of various inclinations. 



3. Inclination Dependency of Disk Composite Spectrum 



A composite spectrum of the disk galaxies is constructed in each of the bins of inclina- 
tion. The flux density (f\) of each disk galaxy spectrum is converted to luminosity density 
(La) by multiplying with the standard factor Amd\, where di is the luminosity distance of 
the galaxy, so that a fair comparison among the composites can be u ndertaken regardless o f 
the redshift of the galaxies. This importance was also pointed out in iBurstein et al.l ( )199ll ). 
We use the concordance cosmological parameter^ in calculating d^. The flux density of 
a composite spectrum in each wavelength bin is taken to be the geometric mean of those 
from the contributing sp ectra, which preserves th e average of the steepness of the spectral 
slopes in a sample (e.g., IVanden Berk et al.l 120011 ). and is effective in removing noisy data. 
Following the SDSS convention, the composite spectra are expressed in vacuum wavelengths, 
and are rebinned into 3700 — 8000 A with a spectral resolution of 70 km s _1 per wavelength 
bin. 

The comparison between the composite spectra of face-on (b/a = 0.9 — 1.0) and edge-on 
(b/a = 0.1 — 0.2) disk galaxies is shown in Figure E] and [71 A marked difference is seen in the 
luminosity and the slope between the continua, as such the inclined galaxies show a lower 
luminosity and a smaller slope in the continuum. The effect is in agreement wit h the usual 
power-law wavelength dependence in extinction curves (e.g.. ICalzetti et al.lll994l ). 



4 The cosmological parameters are IV = 0.73, = 0.27, and h = 0.71. 
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4. Continuum Extinction in g, r, i bands 

4.1. Relative Continuum Extinction 

Assuming the inclination dependency in the composite spectra is mainly due to extinc- 
tion in the samples of the galaxies, we calculate the extinction of the inclined disk at the 
filter band x, A x (b/a) (mag), by comparing the inclined with face-on luminosity (L x (b/a) 
and L x (l)) derived from the composite spectra according to 



A x (b/a) 



AM) 



-2.51og 10 [L x (b/a)/L x (l)] 



(2) 



where A x (l) is the face-on extinction. The corresponding relative optical depth (r x (b/a) — 
t x (1)) can be obtained by dividing the extinction with the factor 2.5 log 10 (e) or 1.086. Al- 
though one can get only the inclined extinction relative to the face-on value but not both 
individually, the advantage of this approach is its being empirical. Throughout the paper, 
we denote the face-on inclinations to be b/a = 1 for simplicity, where the highest inclinations 
are in fact 0.94 ± 0.02 for our sample (Table [2]). 

The broad-band luminosity is derived by convolving a composite spectrum with a given 
filter response curve. To focus on the extinction in the stellar populations instead of other 
components such as the H II regions, the flux contribution from any emission line to the 
synthetic magnitudes are excluded. To make sure the continuum extinction is not biased by 
the continuum estimation method, we use and compare two independent approaches: (1) an 
empirical sliding window approach (2) a Bayesian approach based on a stellar population 
model. 

In the sliding window approach, the continuum luminosity density at a given restframe 
wavelength is estimated to be the average luminosity density between the 40th and 60th 
percentiles of the sample distribution of the contributing luminosity densities, which are in 
turn taken from within ±150 wavelength bins (in the spectral resolution of 70 km s _1 per 
bin) around the centered wavelength. The emission lines are masked out within a constant 



wavele ngth window of ±280 km s 1 , utilizing the line list in Table 30 of IStoughton et al. 
(booj ). 



In the Bayesian approa c h, the method outlined in iKaviraj et al.l (120071 ) is adopted. We 
use the iBruzual fc Charlotl ( 20031 ) high- resolution model (1 A resolution in the restframe 
wavelengths considered in this work) to deter mine the integrated s tellar light of each com- 
posite spectrum, and the extinction model by ICalzetti et al.l (120001 ) for the dust reddening. 
The star formation history is characterized by an exponential declining star formation rate 
defined by an e-folding time, and the model parameters are, the age of the oldest stars: 
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1 — 13.7 Gyr, the stellar metallicity (Z)\ 0.0004 — 0.05, the e-folding time: 1 — 15 Gyr, and 
the nebular gas emission-line reddening: — 0.8 mag. 

For space saving, not all of the results using both continuum estimation methods are 
presented. However, we stress that there is no substantial difference detected in the con- 
tinuum extinction values between both approach. For unity, we present and discuss our 
results in this work based on the Bayesian continuum estimation, unless otherwise specified. 
One significant advantage of the Bayesian approach is in the Balmer-line measurements, as 
will be shown in £j5j where the flux contribution from the underlying stellar absorption is 
represented automatically by using the best-fit stellar population model spectrum. 

In Figure E] we see that the relative g-band extinction A g (b/a) — A g (l) increases with the 
inclination of the galaxies. The relative extinctions for the volume-limited sample are given 
in Table In the highly inclined objects (b/a = 0.26 ± 0.03) A g (b/a) - A g {l) is « 0.4 mag, 
A r (b/a) - A r {l) is « 0.2 mag, and A { (b/a) - A { (1) is « 0.1 mag. 



4.2. Empirical Extinction Model 



To extrapolate the relative extinction to even higher inclin ation values where no result 



is available, empirical models are use d. The log(6/a) model (jde Vaucouleurs et al.l Il991 
Bottinelli et~aDll995l : ITuIIv et al.lll998h is 



A x (b/a) 



AJ1) 



-7* lo Sio 



(3) 



The best-fit model is drawn in Figure [HJ As found also by lUnterborn fc Rydenl (120081 ). the 
actual dependence is found here to be steeper than that described by the log(6/a) model. 
Instead, we found the following 



A x (b/a) -A x (l)=7] x logi (b/a) 



(4) 



provides a better fit to our result. This choice is purely empirical, as a physically motivated 
model depends on knowledge of the actual distribution of stars and dust in a galaxy. In the 
g band, the best-fit rj g =1.2 mag. This value gives a relative extinction of 1.2 ma g for highly 
inclined galaxies with b /a = 0.1, agrees with those from the previous studies (jShao et al. 
20071 : iMaller et al.l 120091 ) at the same inclination, both being > 1 mag in the g band. 
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4.3. Theoretical Extinction Model: Face-On Extinction 

To derive the face-on (b/a = 1) extinction, the inclination-dependent continuum extinc- 
tion is fitted with theoretical models which describe the distributions of stars and dust in 
individual galaxies. The screen model is 



A x (b/a) = A x {l)/{b/a) . 



(5) 



In this scenario, the absorbing layer of material is located above the stars like a screen. 
This choice is motivated by the study of IShao et al.l (120071 ) where the authors found that 
the optical depth of inclined disk galaxies in the SDSS is proportional to the cosine of the 
inclination angle, which is in accord to the description by the screen model. Next, the slab 
model is considered 



A x (b/a) = -2.5 log 



10 



b/a 



-T X {i)/{b/ a y 



(6) 



for which the stars and the dust are mixed uniformly in a slab. Finally, the sandwich model 
is 



A x (b/a) 



-2.51og 10 

C * (b/a 



+ 



\ <L (l + e -r«(l)/(6/«) 



) 



(7) 



(1 



-r x (l)/(b/a) 



for which a layer of dust+stars mixture is sandwiched in-between two layers of stars. In these 
models the face-on extinction and the optical depth are respectively AJ l) and t x (1). The 
geometry of these models are illustrated in detail in iDisney et al.l ( 1989 ). Our model choice 
for dust+stars is non-exhaustive, there are other models in the literature which address, e.g., 



the presence of clumped dust (e.g. . Tuffs et al.ll2004j ). which we plan to explore in the future. 



The best-fit models are drawn in Figure |9j The corresponding best-fit face-on extinctions 
in the SDSS g, r, i bands are given in Table El Among all models the screen model performs 
the least satisfactory, whereas the slab and sandwich models pe rform similarly. Th e best- 
fit values of A g (l) are model-dependent, as pointed out also by IDisney et al.l ( 119891 ). The 
best-fit A g (l) however are similar in both the slab and sandwich models, ~ 0.2 mag. 
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4.4. Extinction Curves in the Optical 

The wavelength dependence of the relative extinction is obtained by comparing the 
inclined to face-on composite spectra. In the following discussions, the inclined spectrum is 
chosen to be that from the highest inclinations (b/a = 0.1 — 0.2), the detailed inclination 
dependence of the relative extinction in the continuum can be obtained, e.g., by invoking the 
best-fit empirical or theoretical models ( §4.2[ §4.3p . Within the restframe wavelength range 
3700 — 8000 A, the relative extinction as a function of wavelength is calculated according 
to Eqn. [2j where the broad-band luminosity on the left-hand-side is now substituted by the 
luminosity density (L\) of the composite spectrum multiplied by the size of the wavelength 
bin (AA) at the concerned wavelength. The resultant extinction curve is shown in Figure [TD] 
for A = 3700—5000 A, and FigurefTTlfor A = 5000—8000 A. Comparing the relative extinction 
in the lines and in the stellar continuum, the following cases can be visually identified: 

f lines which exhibit relative extinction larger than that in the surrounding stellar contin- 
uum: CaK and H absorptions, NaDA5896 absorption, He IA3889, Ha, H/3, H7, H<5 
emissions, [N II]A6550 and A6585 emissions 

f lines which exhibit relative extinction smaller than that in the surrounding stellar contin- 
uum: [O III]A4960 and A5008, [O II]A3727 and A3730 and [O I]A6302 emissions 

I lines which exhibit relative extinction similar to that in the surrounding stellar continuum: 
[S II]A6718 and A6733 emissions 

f uncertain due to small line strength in the composite spectra: Hi], B.6 emissions, [O III]A4364, 
[O I]A6366 emissions 

These lines are all present in the face-on composite disk spectra, albeit in the last case the 
line strengths in the composite spectra are too small for identification. It is intriguing to see 
that there is a variety of behaviors, in terms of the inclination-induced dust extinction among 
the emission/absorption lines, and their relation to the extinction in the stellar continuum. 
In particular, none of the oxygen forbidden emission lines shows larger relative extinction 
than in the stellar continuum. Besides, the relative extinction in the [O III]A5008 is smaller 
than that in the H/3, a little bit surprising if they are emitted from the same H II region. 
These results are expected to give us hints on the geometry of the dust in relation to the 
regions where these lines arise, and on the extinction mechanism (e.g., type of dust) on 
each of the elements. This analysis however is beyond the scope of this paper, and will be 
performed in a future work. In this work we focus only on the Balmer decrement, Ha/H/3, 
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commonly used for deriving the dust reddening in H II region within a emission-line galaxy, 
later in §51 

The wavelength-dependent relative extinction in the stellar continuum is next fitted by 
a polynomial function of the 3rd degree 

3 

A x *Ax(b/a) - A a *aa(1) = a i ^ ' ( 9 ) 

3=0 

where the wavenumber v (i.e., inverse wavelength, v = 1/A) is in the unit of inverse micron, 
/jm" 1 . The best-fit coefficients are given in Table HJ 

Assuming the stellar population and dust properties of the edge-on galaxies are not 
specially different from those of the face-on galaxies, the combination of Eqn. [2]and Eqn. Ocan 
be used to correct against the effects of reddening and extinction in the observed continuum 
of an inclined star-forming disk galaxy. 



5. Balmer Decrement 

The reddening on the emissions from the H II regions such as the Balmer decrement (the 
inten sity ratio between Ha and H /3) depends on the detailed geometry of the dust and the gas 



(e.g. 



Caplan fc Deharvenglll986l ). If the reddening induced by the galaxy inclination affects 
the emissions from the H II regions of the galaxy in a similar manner as the integrated light 
from the continuum-generating stars, the measured Ha/H/3 ratio can in-principle depend 
on the inclination of a galaxy. We therefore explore the relationship between the Balmer 
decrement and the inclination-dependent continuum extinction. 

The Balmer line ratio Ha/H/3 provides a mean for characterizing the dust reddening in 
"normal" H II regions (case B recombination, n e = 100 cm -3 , T e = 10,0 00 K), because th e 



line ratio is fairly in-sensitive to the electron temperature (see Table 4.4 of lOster brocklll989l ) . 
In measuring the EW of a Balmer line, we fit to the continuum-subtracted composite spec- 
trum a single Gaussian function. In measuring the uncertainty in an EW, the following steps 
are performed. Firstly, the uncertainty in the luminosity densities of the stellar continuum in 
each composite is calculated by using 5L C X /L C X = 5L\/L\, where S denotes the uncertainty of 
a quantity. The uncertainty in each Balmer EW is calculated to be J2\ R — L x )/L c x dX, 
where A# is the region of influence of the emission lines, taken to be around the line center, 
±280 km s _1 and ±200 km s _1 , respectively for the emission lines Ha and H/3. Finally, the 
uncertainty in the Ha/H/3 is propagated by using (Ha/H/3) v/((5H/3/H/3) 2 + (5Ha/Ha) 2 ). 
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Examples of the best-fit stellar continua in the vicinity of H/3 and Ha emissions are 
shown respectively in Figure [12] and Figure [13j Using the Bruzual and Chariot stellar 
population model and the Calzetti dust model, in the Bayesian approach the underlying 
stellar absorption lines are represented obviously in the best-fit continuum. 



The inclination dependency of the Ha and H/3 EW's is shown respectively in Figured 
and [151 f° r t ne sliding window approach (left figures) and the Bayesian approach (right 
figures). The offset between both approaches is ~ 7 A for the Ha, and » 4A for the H/3. 
If the stellar absorptions are not taken into account as in the sliding window approach, the 
Ha/H/3 would be mis-taken to be approximately a factor of 1.5 larger for our galaxies. 

The inclination dependency of the measured Ha/H/3 is shown in Figure [TBI using the 
geometric mean spectra. The Balmer decrements are also calculated in th e median spectra, 



show n in Figure [T7] which is expected to preserve the relative line strengths (IVanden Berk et al 



20011 ) . The Ha/H/3 calculated using the median composite spectra are consistent with those 
obtained by using the geometric composite spectra. Both plots show that the Ha/H/3 remains 
fairly constant with inclination - a different behavior from the stellar continuum extinction, 
which shows inclination dependency. 



5.1. Reddening in H II Region vs. in Continuum-Generating Stars 



We adopt the Balmer optical depth (r^, ICalzetti et al.lll994j ). a measure of the dust 
reddening from the measured Balmer line ratio, as follows 



T,, 



Ha/H/3 
(Ha/H/3) T 



(10) 



The superscript / was used by the authors to indicate that derived from the emission lines 
Ha and H/3. The ratio (Ha/H3) T is the un-reddened Balmer decre ment, taken to be the 
theoretical value, 2.86 (jCaplan &: Deharveng||l986l ; lOsterbrockl 119891 ) in this work, which is 
accurate to 1 % for T e = 5, 000 - 20, 000 K. 

The Balmer optical depth vs. the inclination-dependent g-band stellar continuum op- 
tical depth {jg{b/a)) is shown in Figure [TBI The face-on optical depths (r 5 (l)) are taken to 
be the best-fit values in the slab model (Table [3]). One relation is evident, that most of the 
points of small optical depths deviate by more than la sample scatter from the one-to-one 



relationship (the dashed line in Figure [TBJ where r, 
larger. 



Tg(b/a)), because is relatively 
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Next, the color excess obtained from the Balmer decrement is compared with that 
derived from the inclination-dependent continuum extinction. The intrinsic color excess 
derived from the Balmer decrement, labeled as E(B — y ") Balmer in this work, is related to 
the Balmer optical depth as follows ( ICalzetti et al.lll994l ) 



E{B -V) BaW « 0.935 r B . 



The color excess in stars, E(B — V) b J tars , is calculated at a given inclination, b/a, by using 



E(B-Vf^ s = (B-Vf^-iB-Vt^ (12) 
= A B (b/a)-A v (b/a) , (13) 

where (B — V)® taTS is the intrinsic stellar Johnson B — V color of the stars, and (B — V)H® TS 
is the reddened color due to the inclination effect. Looking at Eqn. [13] and Eqn. [2j the value 
of E(B — V) b s (® rs can be calculated by using the empirical continuum extinctions A B (b/a) — 
Ab{1) and Ay(b/a) — Ay(l), obtained in a similar fashion as in §H by using the Johnson 
filter response curves. The face-on extinctions Ab{1) and ^4y(l) are taken to be those from 
the slab model (Table [3]). The E(B — V) stars depends weakly on the model, because the 
model-dependency of the face-on extinctions in the two photometric bands is reduced in the 
subtraction when calculating E(B — V) stars . If the face-on extinction derived by the screen 
model is used instead, the points are to be moved toward left horizontally by 0.1 unit in the 
plot (not shown), which do not affect the above results qualitatively. 

The comparison between the Balmer-line and stellar continuum color excess is plotted 
in Figure dHJ Most of the points deviate by more than la uncertainty from the "uniform 



interstellar extinction" line ( jCaplan fc Deharvenglll986l ). the dashed line in each sub- figure 



which indicates the presence of uniform interstellar extinction in every galaxy of the sample. 
And, the Balmer-line color excess is larger than the continuum color excess. 

If we assume the color excesses are related so that E(B — V)^ B = cE(B — ^) Balmer , 
where c is a constant, then we find 



E{B - V) 



b/a 
stars 



(0.23 ±0.20) E(B-V) 



Balmer 



(14) 



for our sample. The constant is calculated using c = J] E(B — ^Ostars / S E{B — ^) Ba imer' 
over all of t he inclinations. The r esult E[B — ^) Balmer > E[B — V) b s (^ rs agrees qualitatively 
with that by lCalzetti et al.l (120001 ) (their Eqn. 3, with a proportional constant of 0.44±0.03), 



- 16 - 



although our proportional constant is approximately a factor of 2 lower, both values are 
consistent to la uncertainty. We discuss these results in §6.11 



6. Discussion & Summary 

6.1. Spatially Non-uniform Interstellar Extinction 

Our study shows the lack of a correlation between the Balmer optical depth and the 
continuum optical depth, and between the Balmer-line color excess and the continuum color 
excess. If every galaxy of the sample exhibits uniform interstellar extinction, one would 
expect E{B — ^) Balmer = E{B — V) stars at a given inclination. These results are likely the 
manifestations of the non-uniform interstellar extinction in at least some of the galaxies in 
each of the sample. Found also earlier in the H II regions in the Large Magellanic Cloud 



( jCaplan fc Deharvengl 119861 ). it describes the situation where the optical thickness of the 
dust in front of the light-emitting regions is not the same across the area being observed on 
a galaxy. Combining with the factor in which the spatial distribution of the star-forming 
H II regions can be different (e.g., clumped and/or patchy) from that of the continuum- 
generating stars (e.g., low mass stars), the measured extinctions in the stars and in the H II 
regions thus can be different. 



6.2. Summary 

Our study provides a determination of the wavelength-dependent extinction in a large 
and well-defined sample of star-forming disk galaxies, independent from previous works in 
the literature. We found that both the luminosity and the steepness of the spectral slope of 
the composite spectra decrease (i.e., extinguished and reddened) with the inclination of disk 
galaxies, in a volume-limited sample in z = 0.065 — 0.075 and M r = -19.5 to -22 mag from the 
SDSS DR5. Assuming this effect is mainly due to intrinsic extinction in the galaxies, in the 
SDSS g-band the inclined relative to face-on continuum extinctions are found empirically to 
be 1.2 mag for highly inclined objects (b/a = 0.1). The derived face-on extinction is model 
dependent, nonetheless with both the slab and sandwich models giving 0.2 mag in the g band. 
The extinction curve shows a variety of extinction, in both the sign and the amplitude, for 
emission/absorption lines relative to the stellar continuum. The parameterized extinction 
curves allow for correcting the observed galaxy continuum against the inclination effect. 

The Ha /HP does not show dependence with the inclination, and remain fairly constant. 
This inclination dependency is different from that in the stellar continuum, suggesting the 
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mechanism and/or the dust configuration responsible for the reddening are different in the 
H II region and in the continuum-generating stars. The color excess E(B — V) 



B aimer 



> 



E(B — U) stars at most inclinations, rather than the equality of both. From this we conclude 
that at least some of the galaxies exhibit the spatially non-uniform interstellar extinction, 
possibly caused by the presence of the clumped dust in the vicinity of the H II region in the 
galaxies. 

The next s tep is to correct for the observed properties w i th existing parameter-inclination 



relati ons (e.g., IShao et all 120071 : lUnterborn & Rvdenl BOOS Driver et all 120081 : Mailer et al. 



20091 . and this work), re-select the galaxies, and repeat the analysis. The convergence of 
the parameter-inclinat ion relations i s essen tial for ensuring their correctness. This approach 
is with reference with iDriver et all (120071 ). who developed an iterative approach to derive 
cosmic mass densities of the bulge, disk and dust components due to disk galaxies, and found 
converged extinction-inclination relations for those galaxies. 

Our study also raises several interesting questions. For example, does the variation of 
extinction as a function of inclination from our current empirical measurement agree with 
that by other approaches, such as stellar population synthesis modeling? Can we build a 
model to describe the distributions of dust, gas, stars and line-emitting regions, that would 
explain the relative amount of extinction in the emission/absorption lines and the stellar 
continuum, to meet the extinction curve? Can we apply the current analysis to other galaxy 
types? What is the effect of galaxy inclination on photometric redshift, given that the 
inclination causes a reddening to the colors of a galaxy? We will leave these studies to future 
papers. Because our sample selection is well defined in a publicly available data set, works 
can be carried out easily even by other researchers in exploring any secondary factors that 
may affect the extinction in the galaxies. 



7. Acknowledgments 

CWY thanks Reynier F. Peletier for discussions on extinction in galaxies. We thank the 
referee for informative and insightful comments, particularly on the sample selection. We 
acknowledge support through grants from the W.M. Keck Foundation and the Gordon and 
Betty Moore Foundation, to establish a program of data-intensive science at to the Johns 
Hopkins University. IC and LD acknowledge support from NKTH:Polanyi and KCKHA005 
grants. 

This research has made use of data obtained from or software provided by the US 
National Virtual Observatory, which is sponsored by the National Science Foundation. 



-18 - 



Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, 
the Participating Institutions, the National Science Foundation, the U.S. Department of 
Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, 
the Max Planck Society, and the Higher Education Funding Council for England. The SDSS 
Web Site isjhttp://www.sdss.org/ 

The SDSS is managed by the Astrophysical Research Consortium for the Participating 
Institutions. The Participating Institutions are the American Museum of Natural History, 
Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western 
Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Ad- 
vanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute 
for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the 
Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National 
Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for 
Astrophysics (MPA), New Mexico State University, Ohio State University, University of 
Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Obser- 
vatory, and the University of Washington. 

REFERENCES 

Bailin, J., & Harris, W. E. 2008, ApJ, 681, 225 

Beckman, J. E., Peletier, R. F., Knapen, J. H., Corradi, R. L. M., & Gentet, L. J. 1996, 
ApJ, 467, 175 

Boselli, A., & Gavazzi, G. 1994, A&A, 283, 12 

Bottinelli, L., Gouguenheim, L., Paturel, G., & Teerikorpi, P. 1995, A&A, 296, 64 

Burstein, D., Haynes, M. P., & Faber, M. 1991, Nature, 353, 515 

Bruzual, C, & Chariot, S. 2003, MNRAS, 344, 1000 

Calzetti, D., Kinney, A. L., & Storchi-Bergmann, T. 1994, ApJ, 429, 582 

Calzetti, D., Armus, L., Bohlin, R. C, Kinney, A. L., Koornneef, J., & Storchi-Bergmann, 
T. 2000, ApJ, 533, 682 

Caplan, J., & Deharveng, L. 1986, A&A, 155, 297 

Davies, J. I., Phillips, S., Boyce, P. J., & Disney, M. J. 1993, MNRAS, 260, 491 



- 19 - 

Davies, J. L, & Burstein, D. 1995, NATO ASIC Proc. 469: The Opacity of Spiral Disks 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., Jr., Buta, R. J., Paturel, G., & 
Fouque, P. 1991, Volume 1-3, XII, 2069 pp. 7 figs.. Springer- Verlag Berlin Heidelberg 
New York, 

Disney, M., Davies, J., & Phillipps, S. 1989, MNRAS, 239, 939 

Driver, S. P., Popescu, C. C, Tuffs, R. J., Liske, J., Graham, A. W., Allen, P. D., & de 
Propris, R. 2007, MNRAS, 379, 1022 

Driver, S. P., Popescu, C. C, Tuffs, R. J., Graham, A. W., Liske, J., & Baldry, I. 2008, ApJ, 
678, L101 

Elmegreen, D. M. 1980, ApJS, 43, 37 

Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shimasaku, K., & Schneider, D. P. 1996, 
AJ, 111, 1748 

Giovanelli, R., Haynes, M. P., Salzer, J. J., Wegner, G., da Costa, L. N., & Freudling, W. 
1994, AJ, 107, 2036 

Gunn, J. E., et al. 1998, AJ, 116, 3040 

Hogg, D. W., Finkbeiner, D. P., Schlegel, D. J., & Gunn, J. E. 2001, AJ, 122, 2129 

Holmberg, E. 1958, Meddelanden fran Lunds Astronomiska Observatorium Serie II, 136, 1 

Holwerda, B. W., Keel, W. C, & Bolton, A. 2007, AJ, 134, 2385 

Hopkins, A. M., et al. 2003, ApJ, 599, 971 

Huizinga, J. E., & van Albada, T. S. 1992, MNRAS, 254, 677 

Ivezic, Z., et al. 2004, Astronomische Nachrichten, 325, 583 

Kauffmann, G., et al. 2003, MNRAS, 341, 33 

Kaviraj, S., Rey, S.-C, Rich, R. M., Yoon, S.-J., & Yi, S. K. 2007, MNRAS, 381, L74 

Kewley, L. J., Groves, B., Kauffmann, G., & Heckman, T. 2006, MNRAS, 372, 961 

Lupton, R., Gunn, J. E., Ivezic, Z., Knapp, G. R., & Kent, S. 2001, Astronomical Data 
Analysis Software and Systems X, 238, 269 



-20 - 

Mailer, A. H., Berlind, A. A., Blanton, M. R., k Hogg, D. W. 2009, ApJ, 691, 394 

Miller, C. J., Nichol, R. C, Gomez, P. L., Hopkins, A. M., k Bernardi, M. 2003, ApJ, 597, 
142 

Mollenhoff, C, Popescu, C. C, k Tuffs, R. J. 2006, A&A, 456, 941 

Osterbrock, D. E. 1989, Astrophysics of gaseous nebulae and active galactic nuclei, University 
Science Books. 

Pier, J. R., Munn, J. A., Hindsley, R. B., Hennessy, G. S., Kent, S. M., Lupton, R. H., & 
Ivezic, Z. 2003, AJ, 125, 1559 

Ryden, B. S. 2004, ApJ, 601, 214 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 

Shao, Z., Xiao, Q., Shen, S., Mo, H. J., Xia, X., k Deng, Z. 2007, ApJ, 659, 1159 

Shimasaku, K., et al. 2001, AJ, 122, 1238 

Smith, J. A., et al. 2002, AJ, 123, 2121 

Stoughton, C, et al. 2002, AJ, 123, 485 

Strauss, M. A., et al. 2002, AJ, 124, 1810 

Tuffs, R. J., Popescu, C. C, Volk, H. J., Kylafis, N. D., k Dopita, M. A. 2004, A&A, 419, 
821 

Tully, R. B., Pierce, M. J., Huang, J.-S., Saunders, W., Verheijen, M. A. W., k Witchalls, 
P. L. 1998, AJ, 115, 2264 

Unterborn, C. T., k Ryden, B. S. 2008, ApJ, 687, 976 

Valentijn, E. A. 1990, Nature, 346, 153 

Vanden Berk, D. E., et al. 2001, AJ, 122, 549 

White, R. E., Ill, k Keel, W. C. 1992, Nature, 359, 129 

York, D. G., et al. 2000, AJ, 120, 1579 



This preprint was prepared with the A AS IATgX macros v5.2. 




Fig. 1. — Left: the restframe r-band Petrosian absolute magnitude vs. the restframe u - r 
color for the galaxies in the parent sample. The range of fracDeV r is from — 1. Both the 
magnitude and the color are K-corrected. Right: same as the left, but plotted with only the 
galaxies that lie in a fracDeV r range of — 0.1. Compared with the left figure, now the blue 
galaxies are sampled predominantly, with a contamination by the red sequence of ~ 5%. 
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Fig. 2. — The number distribution of the axis ratio {b/a) of the galaxies in our parent 
sample in the range fracDeV r = — 0.1, at three different slices of redshift. The narrowest 
in the redshifts is designed for the purpose of generating volume-limited samples. Without 



additional magnitude cuts, we see that at lower redshifts (Figure 2(a) ) the face-on galaxies are 
relatively under-sampled. Conversely, at higher redshifts (Figure 2(c)) the edge-on galaxies 
are relatively under-sampled. 
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Fig. 3. — The restframe r-band Petrosian absolute magnitude vs. the redshift of the galaxies 
in the parent sample, where fracDeV r ranges from — 0.1. The lower absolute magnitude 
limit corresponds to an apparent magnitude of 17.77 mag. 
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Fig. 4. — The average inclination of the galaxies (parent sample, fracDeV r = — 0.1), plotted 
in the plane of M r vs. redshift. The corresponding number distribution is shown in Figure |3j 
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Fig. 5. — The r-band apparent half-light radius, r e g b ^ a , as a function of inclination of the 
galaxies in the parent sample that lie in a fracDeV r range from — 0.1. The dashed line is 
the best-fit functional form log 10 (r eff 6/,a ) = log 10 (r^ 1 ) — 0.261og 10 (b/a). 
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Fig. 6. — The comparison between the restframe composite spectra of face-on (b/a = 0.9 — 
1.0, black) and edge-on (b/a = 0.1 — 0.2, gray) star-forming disk galaxies, for the volume- 
limited sub-sample. The flux density per unit wavelength is an "observed" quantity in the 
sense that inclination correction is not applied. 
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Fig. 7. — Continuation of Figure [6] to wavelengths 5000 — 8000 A. 
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Fig. 8. — Dependency of the g-band relative extinction in the disk galaxies on their inclina- 
tion, for the volume-limited sample. The best-fit empirical models are the log model (dashed 
line), and the log 4 model (dot-dot-dot-dashed line). 
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Fig. 9. — Dependency of the g-band relative extinction in the disk galaxies on their incli- 
nation, for the volume-limited sample. The best-fit theoretical models are the screen model 
(dotted line), the slab model (dashed line) and the sandwich model (dot-dot-dot-dashed 
line) . 
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Fig. 10. — The extinction curve obtained by comparing the edge-on (b/a = 0.1 — 0.2) with 
face-on (b/a = 0.9 — 1) composite spectra, according to Eqn. [2J The gray line is the best-fit 
3rd degree polynomial function on the emission-line excluded continuum. See Figure [UJ for 
restframe wavelengths longer than 5000 A. 
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Fig. 11. — Continuation of Figure [TO] to wavelengths 5000 — 8000 A. 
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Fig. 12. — The face-on (b/a = 0.9 — 1.0) composite spectrum (black) and the best-fit stellar 
continuum using the Bayesian approach (gray), in the vicinity of the H/3. 




Fig. 13. — The face-on (b/a = 0.9 — 1.0) composite spectrum (black) and the best-fit stellar 
continuum using the Bayesian approach (gray), in the vicinity of the Ha. 
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Fig. 14. — The inclination dependency of the Ha equivalent width. The sliding window 
approach (left) and the Bayesian approach (right) are used respectively in estimating the 
stellar continuum in the composite spectra. The error bar is lcr statistical uncertainty for 
each parameter. 
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Fig. 15. — The inclination dependency of the H/3 equivalent width. The sliding window 
approach (left) and the Bayesian approach (right) are used respectively in estimating the 
continuum in the composite spectra. The error bar is lcr statistical uncertainty for each 
parameter. 
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Fig. 16. — The inclination dependency of the Balmer decrement for the star- forming disk 
galaxies. The Bayesian approach is used in estimating the stellar continuum in the geometric 
composite spectra. The dashed line indicates the theoretical value of normal H II region. No 
obvious trend of Ha/H/3 with inclination is found. 
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Fig. 17. — Same as Figure [16] but using the median composite spectra. The dashed line 
indicates the theoretical value of normal H II region. The Balmer decrement values are 
consistent to la uncertainty with those in Figure [T6j 
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Fig. 18. — The Balmer optical depth vs. the g-band continuum optical depth. The error bar 
in each axis represents la uncertainty in the derived parameter. The dashed line indicates 
the locus if both optical depths are the same. 
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Fig. 19. — The comparison between the color excess derived from the Balmer decrement 
and from the inclination-dependent extinction in the stellar continuum, respectively. The 
error bar in each axis represents la uncertainty in the derived parameter. The dashed line 
indicates the locus if every galaxy in the sample exhibits uniform interstellar extinction, 
where E(B - V% almcr = E(B - V) stars at a given inclination. 
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Table 1: Number of star-forming disk galaxy spectra in th e volume-limited sample. 

bja & number 



0.94 ±0.03 


387 


0.75 ±0.03 


615 


0.65 ±0.03 


625 


0.55 ±0.03 


621 


0.45 ± 0.03 


678 


0.35 ±0.03 


784 


0.25 ±0.03 


684 


0.17 ±0.02 


214 


0.09 ±0.00 


2 



a The average inclination (± la in the sample scatter) of the galaxies. 



Table 2: Relative extinction 3, in disk galaxies as a function of inclination, in SDSS g, r, i 
bands. 



b/a h 


A g (b/a)-A g (iy 


A r (b/a) -A-(l) 


4(V«)-4(i) 


0.94 ± 0.03 


0.00 ±0.18 


0.00 ±0.12 


0.00 ±0.15 


0.75 ±0.03 


-0.04 ±0.18 


-0.04 ±0.12 


-0.04 ±0.15 


0.65 ±0.03 


-0.05 ±0.18 


-0.05 ±0.12 


-0.05 ±0.15 


0.55 ±0.03 


-0.09 ±0.18 


-0.10 ±0.12 


-0.11 ±0.15 


0.45 ±0.03 


-0.05 ±0.18 


-0.09 ±0.12 


-0.11 ±0.15 


0.35 ±0.03 


-0.04 ±0.18 


-0.12 ±0.12 


-0.16 ±0.15 


0.25 ±0.03 


0.15 ±0.18 


0.02 ±0.12 


-0.06 ±0.15 


0.17 ±0.02 


0.41 ±0.19 


0.23 ±0.13 


0.11 ±0.15 



a The restframe continuum relative extinction value, in magnitude. Emission lines are excluded by using the 
Bayesian approach. 

b The average inclination (± la in the sample scatter) of the galaxies. 
c The relative extinction (± la uncertainty) of the galaxies, in magnitude. 
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Table 3: Best-fit face-on extinction a of star-forming disk galaxies in SDSS g, r, % bands. 



model 


A,(i) b 


reduced x 2 


A-(i) 


reduced x 2 


4(1) 


reduced x 2 


screen 


0.06 ±0.03 


0.36 


0.02 ±0.02 


0.67 


-0.01 ±0.02 


0.47 


slab 


0.20 ±0.11 


0.17 


0.09 ±0.06 


0.47 


0.02 ±0.07 


0.45 


sandwich 


0.15 ±0.05 


0.19 


0.14 ±0.14 


0.53 


0.15 ±0.18 


0.45 



a Thc rcstframc continuum extinction value. Emission lines are excluded by using the Bayesian approach. 
fc The best-fit face-on extinction (± la uncertainty) of the galaxies, in magnitude. 



Table 4: Extinction curve 3, in star-forming disk galaxies, parameterized with 3rd degree 
polynomial 13 



a 


ai 


a 2 


a 3 reduced x 2 


-0.554 ±0.014 


0.564 ±0.023 


-0.057 ±0.012 


0.003 ±0.002 0.15 



"Obtained by comparing between the edge-on (b/a = 0.1 — 0.2) and face-on (b/a = 0.9 — 1.0) composite spectra 
in the restframe wavelengths 3700 — 8000 A. Emission lines are excluded by using the Bayesian approach. 
b A\*A\(b/a = 0.1 — 0.2) — ^1a*aa(1) = Sj=o a o where the inverse restframe wavelength v is in the unit 

of \lYOT 1 . 



